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ABSTRACT 

Context. Pulsations in pre-main sequence stars have been discovered several times within the last years. But nearly all of these pul- 
sators are of 6 Scuti-type. y Doradus-type pulsation in young stars has been predicted by theory, but lack observational evidence. 
Aims. We present the investigation of variability caused by rotation and (y Doradus-type) pulsation in two pre-main sequence mem- 
ry^ ' bers of the young open cluster NGC 2264 using high-precision time series photometry from the CoRoT satellite and dedicated high- 

resolution spectroscopy. 

i-C ' Methods. The variability found using the CoRoT data was combined with the fundamental parameters and chemical abundances 

derived from high-resolution spectroscopy, obtained at the Mc Donald Observatory, to discuss the presence of pulsation and rotation 
in the two NGC 2264 cluster members. Time series photometry of NGC 2264 VAS 20 and NGC 2264 VAS 87 was obtained by the 
CoRoT satellite during the dedicated short run SRaOl in March 2008. NGC 2264 VAS 87 was re-observed by CoRoT during the short 
run SRa05 in December 201 1 and January 2012. Frequency analysis was conducted using Period04 and SigSpec. The spectral analysis 
was performed using equivalent widths and spectral synthesis. 

Results. The frequency analysis yielded ten and fourteen intrinsic frequencies for NGC 2264 VAS 20 and NGC 2264 VAS 87, re- 
spectively, in the range from to 1.5d 1 which are attributed to be caused by a combination of rotation and pulsation. The effective 
temperatures were derived to be 6380+150 K for NGC 2264 VAS 20 and 6220± 150 K for NGC 2264 VAS 87. Membership of the two 
stars to the cluster is confirmed independently using X-ray fluxes, radial velocity measurements and proper motions available in the 
literature. The derived Lithium abundances of log n(Li) = 3.34 and 3.54 for NGC 2264 VAS 20 and NGC 2264 VAS 87, respectively, 
are in agreement with the Lithium abundance for other stars in NGC 2264 of similar r cff reported in the literature. 
^\ \ Conclusions. We conclude that the two objects are members of NGC 2264 and therefore are in their pre-main sequence evolutionary 

stage. Their variability is attributed to be caused by rotation and g-mode pulsation rather than rotation only. Assuming that part of 
their variability is caused by pulsation, these two stars might be the first pre-main sequence y Doradus candidates. 

Key words, stars: pre-main sequence - stars: oscillations - stars: individual: NGC 2264 VAS 20, NGC 2264 VAS 87 - techniques: 
CO ' photometric - techniques: spectroscopic 



X 



1 . Introduction from gravitational contraction, have been reported several times 

within the last years (e.g., Ripepi et al. 2006 Zwintz et al. 2009 ), 
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stars (e.g., Rodriguez et al. 2000 Handler 1999). Pulsations in , . , , 

pre-main sequence (PMS) stars, which gain their energy mainly Th ^ classlc f al (Pf»t-) ™nn sequence y Doradus stars have 

spectral types from late A to F and are known to pulsate in high 
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y Doradus pulsators, but no confirmed members of this group 
have been reported yet. 

The search for y Doradus pulsation in PMS stars is com- 
plicated by the likely presence of rotational modulation due to 
spots in such young objects. Candidates for PMS y Doradus pul- 
sation are located in the same region of the Hertzsprung Russell 
(HR) diagram as the hotter members of the group of T Tauri 
stars. The light curves of T Tauri stars are well-known to show 
regular and irregular variations caused by spots on the surface 
and the interaction of the stars with their circumstellar environ- 
ment (e.g., Alencaret al. 2010 ). y Doradus type pulsation occurs 
on the same time scales as variability caused by spots and sur- 
face inhomogeneities, i.e., with periods on the order of up to few 
days. It is therefore hard to unambiguously distinguish between 
the influence of rotation and pulsation in such stars (see, e.g., 
Uytterhoeven et al. 120101 . 

As y Doradus stars have periods on the order of up to a 
few days, they are not easy to observe from ground. Only with 
dedicated multi-site campaigns it is possible to resolve these fre- 
quencies, both spectroscopically and photometrically. The high- 
precision photometry obtained with the CoRoT satellite allowed 
for the first time to study these stars in more detail. 

In this study we present the analysis of CoRoT light curves 
and high-resolution spectroscopy for two PMS members of the 
young open cluster NGC 2264. 

2. NGC 2264: VAS20 and VAS87 

The young open cluster NGC 2264 (a 2 ooo = 6 /! 41'", fcooo = +9° 
53') has been studied frequently in all wavelength regimes in the 
past. It is located in the Monoceros OBI association about 30 pc 
above the galactic plane and has a diameter of ~39 arcminutes. 
Sung et al. (11997) report a cluster distance of 759 + 83pc which 
corresponds to a distance modulus of 9.40 ± 0.25 mag. 

The age of NGC 2264 can only be determined with a rela- 
tively large error as its main sequence consists only of massive 
O and B stars and stars of later spectral types are still in their 
pre-main sequence phase. Additionally, not all cluster members 
might have formed at the same time and star formation is re- 
ported to still be going on in the cluster (e.g., Flaccomio et al. 
1999). Therefore different values for NGC2264's age are re- 
ported in the literature; they lie between 3 and 10 million years 
(e.g., Sung et al. 2004 Sagar et al. 119861 ). If stars are still being 
formed in NGC 2264, some of its members might even be only 
0.15 million years old (Sung et al. |1997l ). 

NGC 2264 is not significantly affected by reddening and 
differential reddening is negligible as well. The mean value of 
E(B -V) = 0.071 ± 0.033 mag reported by Sung et al. ([T997} 
is in good agreement with values found by other authors (e.g., 
Walker fT956l Park et al. 120001 

The first ever discovered PMS pulsators, the two PMS 
8 Scuti stars V588Mon and V589Mon (Breger et al. [19721 . 
are members of NGC 2264. Since then several studies have been 
already devoted to the search for and analysis of its pulsat- 
ing members, e.g., using data from the MOST space telescope 
(Zwintz etal. l20U9l . 

Using the high-precision time series photometry of all acces- 
sible stars in the field of NGC 2264 obtained with the CoRoT 
satellite, we conducted a homogeneous survey for pulsating 
candidates in the magnitude range between 10 and 16 mag 
in R. Cluster members showing variability on time scales be- 
tween about 0.3 to a few days have been the prime candidates 
to search for PMS y Doradus pulsation. Two of these are 
NGC 2264 VAS 20 (CoID0500007018) and NGC2264VAS87 



(CoID0223979759), named VAS 20 and VAS 87 hereafter, and 
are subject of this work. 

VAS 20 (RA 2 ooo = 06:40:05.66, DE 2 oqq = +09 :35:49.3) has 
a V magnitude of 11.248 mag (Sung et al. 1997 ) and a photo- 
metrically determined spectral type F5 (Young 1978 1. An X-ray 
flux was detected by ROS AT/HRI observations (Flaccomio et al. 
2000) favoring the star's youth. Additionally, Furesz et al. (2006 ) 
listed the star as a cluster member according to its radial veloc- 
ity. The membership probability is given as 94% by Vasilevskis 
et al. (11965b . King (1998) reports an effective temperature (T e ff) 
of 6219 ± 140K using B - V and V - I photometry of Sung 
et al. ( 1 19971 ) and calibrations of Bessel (II 9791 ), and a Lithium 
abundance log «(Li) = 3.56 derived from low-resolution spec- 
troscopy. 

VAS 87 (RA2000 = 06:40:56.95, DE 2 q oo = + 09:48:40.7) has a 
V magnitude of 12.323 mag (Sung et al. 119971 ) and a photomet- 
rically determined spectral type F9 (Young 1978 ). The member- 
ship probability is given as 96% by Vasilevskis et al. ( I19651 ). 
VAS 87 is also an X-ray source detected by the CHANDRA 
satellite (Ramirez et al. 2004). Therefore the star is very likely 
to be young and a member of NGC 2264. In the study by King 
([19981 a r eff = 6135 + 140K and a log n(Li) = 3.65 are given. 

3. The CoRoT light curves 

3.1. Observations and data reduction 

The CoRoT satellite (Baghn l2"0"06T) was launched on 2006, Dec 
27th, from Baikonur aboard a Soyuz rocket into a polar, inertial 
circular orbit at an altitude of 896 km. CoRoT carries a 27-cm 
telescope and can observe stars inside two cones of 10° radius, 
one at RA = 06 : 50 and the other at RA = 18 : 50. The field of 
view of the telescope is almost circular with a diameter of 3.8 °. 
Each CCD is a square of 1.3 °. The filter bandwidth ranges from 
370 to lOOOnm. 

CoRoT observed NGC 2264 including VAS 20 and VAS 87 
for the first time for 23.4 days in March 2008 during the Short 
Run SRaOl within the framework of the Additional Programme 
(Weiss 2006). A second short run, SRa05, on the cluster 
NGC 2264 was conducted in December 2011 / January 2012 
during which only VAS 87 was re-observed with a time base of 
about 39 days. 

For both observing runs, the complete cluster was placed in 
one Exofield CCD and data were taken for all stars in the ac- 
cessible magnitude range, i.e., from 10 to 16 mag in R. The 100 
brightest stars in the field of NGC 2264 were the primary targets 
to search for stellar pulsations among PMS cluster members. 

The reduced N2 data for both stars and both runs were ex- 
tracted from the CoRoT data archive. The CoRoT data reduction 
pipeline (Auvergne et al. 2009) flags those data points that were 
obtained during passages of the satellite over the South Atlantic 
Anomaly (SAA) and replaces them with linearly interpolated 
values. We did not use these 'SAA-flagged' data points in our 
analysis. The light curves of VAS 20 and VAS 87 obtained in 
2008 (see Figures Q] and top left panel in Figure [2} consist of 
2704 and 2600 data points, respectively, with a sampling time of 
512 sec. The 201 1/12 data set of VAS 87 has been observed with 
a sampling time of 32 sec and therefore comprises 90073 data 
points. We checked for the presence of contaminant stars in the 
point-spread function and mask used for the observations for all 
three light curves and found no contaminating sources. 

For VAS 20 chromatic light curves in red, green and blue 
were obtained by CoRoT, while for VAS 87 data were taken only 
in white light. The chromatic light curves were used to check 
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the variability of VAS 20. The variations are present in all three 
colors with the expected different photometric amplitudes. The 
top panels of Figures[T]and[2]show the light curves in white light 
for the two stars where the fluxes obtained in the three colors for 
VAS 20 were summed to increase the signal-to-noise. 



-1 5.680 



VAS 20 was also observed by the MOST (Walker et al. 2003 1 
satellite in 2006 (Zwintz et al. 12009b and VAS 87 was included 
in the MOST observations of NGC 2264 in 201 1/12. As the two 
stars are rather faint targets for MOST, the low amplitude vari- 
ability is completely hidden in the noise and the respective light 
curves cannot be used for any further investigation. 



3.2. Comparison of the CoRoT data of SRa01 and SRa05 

The light curves of VAS 87 (Figure |2]i observed in 2008 (i.e., 
SRaOl) and 201 1/12 (i.e., SRa05) show a dramatic difference in 
the amplitudes of the observed variability. A detailed compara- 
tive study of the characteristics of the NGC 2264 observations by 
CoRoT obtained in both years showed that several instrumental 
effects are present. These instrumental effects, as described be- 
low, make it very difficult for us to interpret amplitude changes 
between the two observing years. For that reason we do not con- 
sider amplitude variations in this paper. 

Different instrumental settings were used for CoRoT in 2008 
and 2011/12 which might affect the measured amplitudes: For 
the observations in 2008, the CCD El was used, but as it is 
not operative any more, the observations in 2011/12 were con- 
ducted using the second Exoplanet CCD, E2. For the short 
run SRaOl the best possible observing season for NGC 2264 
was chosen (i.e., March 2008) minimizing the influence of 
scattered light. Due to the coordinated simultaneous observa- 
tions of CoRoT with the satellit es Spitzer (Werner et al. 120041), 
Chandra (Weisskopf et al. 120021 and MOST (Walker et al. l2003l 
in 2011/12 not the optimum observing period was taken (i.e., 
December 201 1 and January 2012). Therefore different orienta- 
tions of the CoRoT satellite in the two runs had to be chosen and 
different photometric templates had to be used. The different ori- 
entation of the cluster on the CCD could have some impact on 
the properties of the observations, e.g., the CCD temperature in 
the 201 1/12 observations was higher than in 2008. Additionally, 
the aging of the CCD resulted in an increased rate of hot pixels 
and a higher dark current level. 

Many of these effects are corrected by the CoRoT pipe-line 
in a first approximation, and do not affect the time-domain prop- 
erties. But the situation is more difficult concerning the ampli- 
tudes of individual modes. For instance, aperture photometry can 
be affected by the background levels of the stellar field and of the 
CCD. An uncertainty on the estimation of the background level 
of 10% leads to an uncertainty on the amplitude of the same 
order. The estimate of the background used by the standard cor- 
rection in the CoRoT pipe-line depends on many factors, such 
as the proximity of the background windows, the orientation of 
the CCD, and the shape of the photometric template. So it can 
vary from one target to another. A detailed study of these effects 
is currently in preparation (Baglin et al. 120 1 2) . 

As the above described effects influence the ultra-precise 
CoRoT data significantly, the relative amplitudes for VAS 87 ob- 
tained from the SRaOl data are different to those derived from 
the SRa05 data. This is why in this paper we do not discuss am- 
plitude variations. 
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Fig. 1. CoRoT light curve (top panel) for VAS 20, original am- 
plitude spectrum (grey) from to 3 d -1 (bottom X axis), i.e., 
from to 34.72 /uliz (top X axis), with the first three intrinsic 
frequencies identified (black), residual amplitude spectra after 
prewhitening the three highest amplitude frequencies (where the 
remaining 7 intrinsic frequencies are again marked in black). 
Note the different y-axes scales. 



3.3. Frequency Analysis 

For the frequency analyses, we used the software package 
Period04 (Lenz & Breger 2005) that combines Fourier and least- 
squares algorithms. Frequencies were then prewhitened and con- 
sidered to be significant if their amplitudes exceeded four times 
the local noise level in the amplitude spectrum (Breger et al. 
1993 Kuschnig et al. [19971 ). We verified the analysis using 
the SigSpec software (Reegen [20071 1. SigSpec computes signif- 
icance levels for amplitude spectra of time series with arbitrary 
time sampling. The probability density function of a given am- 
plitude level is solved analytically and the solution includes de- 
pendences on the frequency and phase of the signal. 

In 2008, VAS 20 and VAS 87 were observed with a sampling 
time of 512 s, corresponding to a Nyquist frequency of 84 d~'; 
hence, the frequency analysis was conducted in the range from 
to 84 d _1 . The Rayleigh frequency resolution l/T corresponds to 
0.04d~', where T is the length of the observations (23.4days). 
In the case of two frequencies present within l/T, Kallinger et 
al. (f2008 ) have shown using numerical simulations that the up- 
per limit for the frequency error is 1 /(4 • T) (see also their Figure 
3) if the SigSpec significance (sig) is higher than 16. In our case 
this corresponds to a resolution of O.Old . For sig < 16, the 
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Fig. 2. CoRoT light curves (top panels) for VAS 87 from SRaOl (2008 data, left panel) and SRa05 (201 1/12 data, right panel; light 
curve with 32sec sampling in grey and binned to 5 12sec in black) to the same X-axes scales, original amplitude spectra (grey) from 
to 2 d _1 (bottom X axis) i.e., from to 23. 15 fjHz (top X axis), with the six common intrinsic frequencies identified in black (middle 
panels), residual amplitude spectra after prewhitening the six frequencies (bottom panels) where the additional eight frequencies 
found in the 201 1/12 data are marked in black (bottom right panel). Note the different Y-axes scales in the top and middle panels. 



frequency error is 1 /(T * -\fsig), i.e., the values given in paren- 
theses in Tables Q] and [2] In our case only F10 in VAS 20 has 
sig < 16. 

In 201 1, VAS 87 was observed with a sampling time of 32 s, 
corresponding to a Nyquist frequency of 1218 d _1 . The complete 
run lasted for about 39 days, but due to instrumental effects, the 
first few days had to be omitted from the analysis. The total 
length of the data set (shown in Figure [2} used in this inves- 
tigation comprises ~35 days, which corresponds to a Rayleigh 
frequency resolution of 0.028 d -1 . Applying the Kallinger et al. 
(2008) criterion, we get an upper limit of the frequency error 
of 0.007 d -1 , if the SigSpec significance is higher than 16. The 
frequency errors for sig < 16 are given in Table [2] 

Using the CoRoT light curve from 2008, VAS 20 shows ten 
intrinsic frequencies in the range from to 1.35 d -1 (Figure Q] 
and TableQ}. Linear combinations between the frequencies have 
been sought within three times the frequency uncertainty. Three 
frequencies, i.e., F4, F9 and F10, might be explained as linear 
combinations, the other seven are clearly independent from each 
other. We believe that the frequency corresponding to the rota- 
tional period of VAS 20 is Fl (see Section 5). 

For VAS 87 the data sets from 2008 and 201 1 were analyzed 
independently and the results were then compared to each other. 
As the data set obtained in 201 1 is significantly longer, the noise 
level in the amplitude spectrum is decreased and additional fre- 
quencies could be identified. Six frequencies are common to 
both data sets (see Table 0, eight additional frequencies were 



found in the 2011 data set. From these 14 frequencies in total, 
seven could be explained as linear combinations (again within 
three times the frequency uncertainties) which are identified in 
the last column of Table [2] and seven are clearly independent 
from each other. Due to the instrumental differences in the ob- 
servations from 2008 and 2011 (see Section l3~2t . the respective 
amplitudes of the common frequencies cannot be used for an as- 
trophysical interpretation. For VAS 87 we believe that F5 is the 
frequency which indicates the rotational period (see Section 5). 

Applying the criteria by Kallinger et al. (2008 ), all frequen- 
cies for VAS 20 and VAS 87 are resolved and are therefore con- 
sidered intrinsic to the stars. 



4. Spectroscopic analysis 

Fundamental parameters and abundances of both stars were de- 
termined from spectra obtained on 2010, December 15th and 
16th with the Robert G. Tull Coude Spectrograph (TS) on the 
2.7-m telescope of Mc Donald Observatory. In the adopted con- 
figuration the cross-dispersed echelle spectrograph has a re- 
solving power of 60 000 and delivered spectra of VAS 20 and 
VAS 87 with signal-to-noise ratios (S/N) per pixel, calculated at 
~5000A, of 100 and 65, respectively. Each spectrum covered 
the wavelength range from 3633-10849 A with gaps between 
the echelle orders at wavelengths larger than 5880 A. 



4 



K. Zwintz et al.: PMS y Doradus candidates in NGC 2264 



Table 1. Results of the frequency analysis of VAS 20: frequencies and periods, amplitudes, signal-to-noise values and SigSpec 
significances are given sorted by the prewhitening sequence. The respective last-digit errors of the frequencies computed according 
to Kallinger et al. (2008) are given in parentheses. Possible linear combinations are listed in the last column. 



No. 


freq 


freq 


period 


amp 


S/N 


sig 


combi 


# 


[d- 1 ] 


L"Hz] 


[d] 


[mmag] 








Fl 


0.672(2) 


7.78(2) 


1.488 


2.726 


14.8 


428.9 




F2 


0.226(3) 


2.61(4) 


4.427 


0.818 


8.6 


164.6 




F3 


0.117(4) 


1.35(5) 


8.566 


0.638 


7.4 


118.6 




F4 


0.347(4) 


4.02(4) 


2.879 


0.572 


8.9 


136.3 


3F3 


F5 


0.160(5) 


1.85(6) 


6.257 


0.380 


6.5 


72.8 




F6 


0.760(6) 


8.80(7) 


1.315 


0.299 


4.7 


47.4 




F7 


0.397(8) 


4.60(9) 


2.518 


0.217 


4.3 


27.2 




F8 


0.723(9) 


8.37(10) 


1.383 


0.320 


7.5 


22.3 




F9 


0.506(9) 


5.86(11) 


1.977 


0.182 


4.5 


20.3 


F1-F5 


F10 


1.350(13) 


15.62(15) 


0.741 


0.105 


3.6 


11.5 


2F1 



Table 2. Results of the frequency analysis of VAS 87 using data from 2008 and 201 1: frequencies and periods, amplitudes, signal- 
to-noise values and SigSpec significances for both years are given sorted first by the prewhitening sequence of the 2008 data and 
then of the 201 1/12 data. The respective last-digit errors of the frequencies computed according to Kallinger et al. (2008 ) are given 
in parentheses. In the last column possible linear combinations are listed. 



No. 


freq08 


freq08 


amp08 


S/N-08 


sig08 


freq 1 1 


freql 1 


ampl 1 


S/N- 11 


sig 11 


period 


combi 


# 


[d- 1 ] 


L"Hz] 


[mmag] 






[d- 1 ] 


I^Hz] 


[mmag] 






[d] 




Fl 


0.519(3) 


6.01(3) 


1.081 


9.6 


219.3 


0.511(2) 


5.92(2) 


0.451 


5.1 


323.8 


1.925 




F2 


0.309(3) 


3.58(4) 


0.410 


10.6 


167.6 


0.291(3) 


3.36(3) 


0.769 


13.8 


100.1 


3.236 




F3 


0.236(4) 


2.73(5) 


0.641 


9.0 


114.1 


0.251(1) 


2.90(1) 


1.996 


19.9 


748.1 


4.235 




F4 


0.554(7) 


6.41(9) 


0.250 


6.1 


31.9 


0.540(2) 


6.25(2) 


0.404 


6.2 


243.6 


1.806 


F2+F3 


F5 


0.274(9) 


3.17(10) 


0.447 


9.6 


23.1 


0.2620(3) 


3.033(4) 


4.098 


23.6 


7980.9 


3.655 


F1-F3 


F6 
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73.0 
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0.243 
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7.47(5) 
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Bias and Flat Field frames were obtained at the beginning of 
each night, while several Th-Ar comparison lamp spectra were 
obtained each night for wavelength calibration purposes. The re- 
duction was performed using the Image Reduction and Analysis 
Facility^ (IRAF). The spectra were normalized by fitting a low 
order polynomial to carefully selected continuum points. 

To compute model atmospheres of the target stars we em- 
ployed the LLmodels stellar model atmosphere code by Shulyak 
et al (2004), while we calculated synthetic spectra with SYNTH3 
(Kochukhov 120071) . Our main source of atomic parameters for 
spectral lines is the VALD database (Kupka et al. Q999). The 
LTE abundance analysis was based on equivalent widths, ana- 
lyzed with a modified version (Tsymbal 1996) of the WIDTH9 
code (Kurucz [1993l ). 

We determined the microturbulence velocity (i/mk) impos- 
ing the equilibrium between line abundance and equivalent 
width for Fei lines, logg imposing the ionisation equilibrium 
for Fe, and r e ff imposing the excitation equilibrium for Fe i 



1 IRAF (http://iraf.noao.edul is distributed by the National Optical 
Astronomy Observatory, which is operated by the Association of 
Universities for Research in Astronomy (AURA) under cooperative 
agreement with the National Science Foundation. 



lines. For VAS 87 we obtained r eff =6220±150K, logg=3.8±0.2 
and y m i c =1.7±0.5kms _1 , while for VAS 20 we obtained 
r eff =6380+150K, logg=4.0±0.2and ^=1.8+0.5 kms" 1 . The 
excitation and ionisation equilibria obtained for VAS 20 and 
VAS 87 are shown in Figures [3] and |4] For both stars, the mea- 
sured equivalent widths span between about 20 and 250 mA. 
Given the fact that the uncertainty on T e ff is only statistical and 
does not take into account the possible presence of various sys- 
tematic effects, we increased the uncertainty to 200 K. By im- 
posing the ionisation equilibrium, we took into account NLTE 
effects for Fe, as given by Mashonkina et al. (1201 11 1. 

It is extremely important to check the parameters obtained in 
this way with other indicators as this method might be strongly 
affected by systematics (see for example Ryabchikova et al. 2009 
and Fossati et al. 1201 11 1. For this reason we used hydrogen lines 
and lines with extended wings to check and log g, respec- 
tively. Note that at this temperaure the hydrogen lines have the 
advantage of not being very sensitive to gravity variations. 

In PMS stars various spectral lines, particularly the hydro- 
gen lines, could show emission as young stars might still be sur- 
rounded by the remnants of their birth clouds. If the hydrogen 
lines would be partially or completely filled, the determination of 



5 



K. Zwintz et al.: PMS y Doradus candidates in NGC 2264 



og(Fe1/N„) »b. x „,X = 




6300 6400 



|MLTE co„ 



o — 

5 hr- 



Fig. 3. Top panel: correlation between line abundance and exci- 
tation potential for Fe i as a function of T e ff for VAS 20. Middle 
panel: Fe i (asterisks) and Fe n (diamonds) abundance as a func- 
tion of log g for VAS 20. A small horizontal shift was applied to 
the points for visualisation purposes. Bottom panel: excitation 
equilibrium for both Fe i (asterisks) and Fe n (diamonds) lines 
for VAS 20. The correlations are shown respectively by a solid 
and a dashed line. 
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Fig. 4. Same as Figure[3] but for VAS 87. 



the effective temperature (T e s) from these lines would be impos- 
sible. For both VAS 20 and VAS 87 there is no sign of emission 
and if any it would affect just the line cores which cannot anyway 
be taken into account due to NLTE effects. We therefore used the 
Hy line to estimate T e g, Ha was not covered and H/3 is affected 
by a defect in the imaging system of the spectrograph. We nor- 
malized the Hy line continuum using the artificial flat-fielding 
technique described in Barklem et al. (2002J, which has proven 
to be successful with TS data (e.g., Fossati et al. 1201 11 1. To re- 
duce systematic uncertainties in the temperature determination 
due to the hydrogen lines normalization, we obtained low res- 
olution {R -10000) spectra, covering Ha (6470-6710 A), with 
the 1.8-m telescope of the Dominion Astrophysical Observatory 
(Canada). Low resolution spectra allow a better control of the 
normalization. 



Within the uncertainties, the hydrogen lines confirmed the 
effective temperatures we derived for both stars with the exci- 
tation equilibrium. We confirmed the derived log g values with 
the analysis of the Mgi line profile at A 5 172 A (for more de- 
tails on the adopted procedure see Fossati et al. 1201 1) , which, 
in this temperature regime, is particularly sensitive to gravity 
variations. Our fundamental parameters provide the best over- 
all description of the available observables: hydrogen and metal- 
lic line profiles. Our analysis shows that the two stars are very 
similar and this is demonstrated also in Figure [5] which shows a 
comparison between the spectra of the two stars. 




5200 5210 
wavelength [A] 




4340 4342 4341 4346 4348 
wavelength [A] 

Fig. 5. Comparison between the observed spectra of VAS 20 
(dashed line) and VAS 87 (solid line) in the region around the 
Mg ib triplet (top panel) and Hy line (bottom panel). Both panels 
present for comparison a synthetic spectrum (dash-dotted line) 
calculated with r e ff=6300K and logg=3.9. To better appreciate 
the similarities between the two stars, we convolved the spec- 
trum of VAS 87 to the same v sin i of VAS 20. 



By fitting synthetic spectra to several weakly blended 
lines we measured a usin/ of 42±2kms~ 1 for VAS 20, and 
18±lkms~'for VAS 87. Given the relatively large vsmi, we 
also performed the analysis of VAS 20 using spectral synthesis 
(for more details see Fossati et al. 2007 1, obtaining results con- 
sistent with the equivalent width results. Given the T e s of the 
two stars we expect the presence of some macroturbulence ve- 
locity (Vmac)- F° r VAS 20, v sin i is so large that the line profiles 
do not require any Umao while for VAS 87, to be able to fit si- 
multaneously the core and the wings of strong lines, we need a 
fmac of about 12kms~'. This value is slightly larger than what 
is expected for main sequence stars of this r e ff ( Valenti & Fisher 
120051) . which might be connected to the fact that this star is in its 
PMS phase. For VAS 20 we were not able to measure i/ mac be- 
cause of the rather large v sin i which prevented us from seeing 
the effects of i/ mac broadening on the line wings. In this respect, 
the v sin i we obtained for VAS 20 is likely to contain also the 
hidden information on w mac , which cannot be directly measured. 

For each star we determined the abundances of about 15 el- 
ements, which are given in Table [3] Notice that the spectrum 
of VAS 87 is badly affected by spikes, making both the con- 
tinuum normalisation and the equivalent width measurements 
rather uncertain, in particular in the blue spectral region. This 
is responsible for the slightly larger abundance uncertainties for 
VAS 87, compared to that of VAS 20. For both stars, we ob- 
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Table 3. LTE atmospheric abundances in program stars with the 
error estimates based on the internal scattering from the num- 
ber of analysed lines, n. For comparison purpose, the last col- 
umn gives the abundances of the solar atmosphere calculated by 
Asplund et al. d2009l >. 



Ion 


VAS 20 




VAS 87 




Sun 




lOgW/Mnt) 


n 


\os(N/N, nt ) 


n 


loe(AVMnt) 

'"bv * / ■* * tot / 


Li i 


-8.55+0.16 


1 


-8.60+0.10 


1 


-10.99 


Nai 


-5.83±0.08 


2 


-5.68: 


2 


-5.87 


Ms i 


-4.63±0.03 


4 


-4.55+0.12 


4 


-4.44 


Si i 


-4.58±0.09 


7 


-4.87+0.14 


5 


-4.53 


Sin 


-4.49: 


1 






-4.53 


Cai 


-5.64±0.09 


16 


-5.71+0.14 


10 


-5.70 


Sen 


-8.97±0.09 


6 


-9.16: 


1 


-8.89 


Til 


-7.00±0.08 


6 


-7.08+0.15 


3 


-7.09 


Tin 


-7.03±0.06 


8 


-7.14+0.13 


6 


-7.09 


Cri 


-6.37+0.17 


15 


-6.39+0.10 


6 


-6.40 


Crn 


-6.44+0.14 


8 


-6.06+0.08 


5 


-6.40 


Mm 


-6.62+0.15 


4 


-6.44+0.14 


5 


-6.61 


Fei 


-4.55+0.10 


127 


-4.63+0.15 


45 


-4.54 


Fen 


-4.52+0.08 


17 


-4.58+0.20 


14 


-4.54 


Co i 






-7.13: 


1 


-7.05 


Nil 


-5.77+0.10 


16 


-5.83+0.08 


8 


-5.82 


Cui 


-8.03: 


1 


-7.68: 


1 


-7.85 


Yn 


-9.94+0.08 


2 


-9.66+0.03 


2 


-9.83 


Ban 


-9.29+0.04 


3 


-9.51+0.05 


2 


-9.86 



tained abundance values comparable to the solar ones, except 
for Ba and Li which are overabundant, in a reasonable agreement 
within the errors to the overall cluster metallicity of -0. 1 5 (Lynga 
119871 1. For Li we derived an abundance of log n(Li)=3.49+0.16 
for VAS 20 and log n(Li)=3.44+0.10 for VAS 87 (see FigureHJ. 
The difference in the uncertainties of the Li abundance for the 
two stars reflects the difference in the projected rotational ve- 
locity (usin/). We obtained the Li abundance from the lines 
at A ~6707A adopting hyperfine structure from Smith et al. 
( 1998 ) and the meteoritic/terrestrial isotopic ratio Li 6 /Li 7 = 0.08 
(Rosman & Taylor 1998 ). The given errors on the Li abundance 
were calculated taking into account the uncertainties on the fun- 
damental parameters (see Fossati et al. 2009). 

4.1. Spectral energy distribution and Hertzsprung-Russell 
diagram 

Figure [7] and [8] show the fits of the synthetic fluxes, calcu- 
lated with the fundamental parameters derived for VAS 20 and 
VAS 87, respectively, to the observed Johnson (Sung et al. |1997l ), 
2MA SS (Zacharias et al. I2005T), Sp itzer (IRAC - Sung et al. 
I20091 and WISE (Cutri et al. l20T2l ) photometry, converted to 
physical units. For this we used the calibrations provided by 
Bessel et al. ( I1998I ), van der Bliek et al. (I1996I ), Sung et al. 
(2009), and Wright et al. (2010), respectively. Adopting the clus- 
ter distance and reddening given by Sung et al. (1997 ), we esti- 
mated the stellar radius of VAS 20 to be 3.5+0.5 R/Rq, while for 
VAS 87 we derived a radius of 2.2±O.5R/R . This comparison 
also confirms the temperature we obtained for the two stars. 

Along the line of sight of VAS 87 there is a dense cloud of 
material which we believe being behind the star, as it would 
otherwise be clearly visible in the spectral energy distribution, 
which we can fit with the same reddening applied for VAS 20. 
This is also supported by the presence of only small amounts 




6704 6706 6708 6710 6712 

wavelength (A) 



Fig. 6. Comparison between the observed spectra (black solid 
line) of VAS 20 (upper panel) and VAS 87 (lower panel) in the 
region of Li line at A ~6707 A and synthetic spectra (red dashed 
line) calculated with the best fitting Li abundance. Note the dif- 
ferent y-axes scales. 




10* 10 s 
wavelength [A] 

Fig. 7. Comparison between LLmodels theoretical fluxes (full 
line), calculated with the fundamental parameters derived for 
VAS 20, with Johnson (crosses), 2MASS (diamonds), Spitzer 
(triangles), and WISE (squares) photometry converted to physi- 
cal units. 



of reddening across the cluster reported in the literature (e.g., 
Sung et al. 1 19971 ). Figure[8]shows also a steep increase in flux in 
the WISE photometry of the two red channels, but this feature is 
not present in the Spitzer photometry. It is therefore possible that 
the WISE photometry has been contaminated by the background 
cloud which would be rather bright at infrared wavelengths. 

Using the spectroscopic temperatures, we put the two stars in 
the HR-diagram. We calculated the luminosity of each star pho- 
tometrically, adopting the magnitudes in the V band, a cluster 
distance of 759+8 3 pc and a reddening of E(B- V)=0.07 1 +0.033 
(Sung et al. 119971 ). We used the bolometric correction by Balona 
([19941 For VAS 20 and VAS 87 we obtained log L/L = 1.235 
and 0.805, respectively. Given all the uncertainties and the fact 
that we adopted a bolometric correction valid for main-sequence 
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wavelength [A] 

Fig. 8. Same as Figure El but for VAS 87. 

stars, we decided to apply a conservative uncertainty on the lu- 
minosity of 0. 1 dex. 
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Fig. 9. Positions of VAS 20 (asterisk) and VAS 87 (diamond) in 
the HR-diagram in comparison with evolutionary tracks (bot- 
tom panel) and isochrones (top panel). In both panels the ZAMS 
is indicated by a dotted line. Evolutionary tracks for PMS 
stars with masses of 2.25 M/M , 2.00 M/M , 1.75 M/M , and 
1.5OM/M are shown respectively by a solid, dashed, dash- 
dotted and dash-triple dotted line. Isochrones with ages of 
lOMyr, 8Myr, 6Myr, and 2Myr are shown respectively by a 
solid, dashed, dash-dotted and dash-triple-dotted line. The bor- 
ders of the y Doradus instability strip (Dupret et al. 2004} are 
marked with thick solid lines. 

Figure [9] shows the position of the two stars in the HR- 
diagram, in comparison to PMS evolutionary tracks from 
Guenther (2012 priv. communication) using the YREC evolu- 
tion code (Demarque et al. 2008) with physics as described in 
Guenther et al. (2009) and isochrones by Demarque et al. (2008). 
For both evolutionary tracks and isochrones we adopted solar 
metallicity. Also included in Figure [9] is the instability strip for 
y Doradus stars (Dupret et al. 120041 . 

The derived values are relatively low compared to the 
theoretical computations of the instability strip for PMS y 



Doradus stars (see Bouabid et al. 1201 11 their Figure 8). The 
cool border of the predicted theoretical (PMS) y Doradus in- 
stability strip lies at about 6700 K, but such a high T e s would 
not allow to fit both hydrogen and metallic lines for both stars. 
Therefore the two stars lie slightly outside the predicted, theoret- 
ical PMS y Doradus instability strip (Bouabid et al. 1201 It . This 
is not too surprising as similar discrepancies are found in other 
groups of (post-)main sequence pulsators. In particular several 
of the y Doradus stars found with data from the Kepler satellite 
(Koch et al. 120101) are located at temperatures similar to that of 
VAS 20 and VAS 87 (Balona et al. l20T2l . 

Additionally, the theoretical seismic properties of PMS y 
Doradus have been investigated without a confirmed member of 
this group of pulsators. A comparison of the observations to the 
PMS y Doradus theory might yield new constraints on the lo- 
cation of the respective instability region. More theoretical and 
observational work will be performed in the future to investi- 
gate this in detail. As the PMS y Doradus instability strip coin- 
cides well with the instability strip of the (post-) main sequence 
y Doradus stars, in Figure [9] we show the location of VAS 20 
and VAS 87 in the HR-diagram with respect to the "classical" 
y Doradus instability strip by Dupret et al. (2004). 

On the basis of the evolutionary tracks we derived masses 
and radii for the two stars. For VAS 20 we obtained a mass 
of 2. 17+0. 15 M , a radius of 3.5+0.5 R Q and a surface grav- 
ity of 3.69+0.11, while for VAS 87 we obtained a mass of 
1.55+0. 15 M , a radius of 2.1 ±0.5 R Q and a surface gravity of 
3.98+0. 1 1 . The stellar radii, derived from the evolutionary tracks 
are in remarkable good agreement with what we obtained from 
the fit of the spectral energy distribution. Within one-sigma, the 
spectroscopic log g values are also in agreement with those de- 
rived from the evolutionary tracks. 

VAS 20 and VAS 87 are two stars born from the same molec- 
ular cloud that are members of NGC 2264 and that have nearly 
the same T e g and log g values, but slightly different masses, radii 
and luminosities (see Table|4]for an overview of the fundamental 
parameters) which place them at different locations in the HR- 
diagram. At present we can only speculate about the reason for 
this discrepancy. 

One possible explanation might be connected to the pres- 
ence of sequential star formation in NGC 2264 and that its star 
forming activity is still going on (e.g., Flaccomio et al. 1999). 
The ages of stars in the cluster lie between 0.8 and 8 million 
years where intermediate mass stars and very low mass stars 
were formed first and and massive stars down to about one so- 
lar mass were formed later on (Sung et al. 11997) . To assess if 
an age spread in NGC 2264 would be a possible reason for the 
slight differences in the fundamental parameters of the two stars, 
we included the isochrones in Figure [9] and obtained an age of 
4±1 Myr for VAS 20 and 10+2 Myr for VAS 87. In the present 
case, the higher mass star (VAS 20) would have to be formed 
later than the lower mass star (VAS 87). This would fit to the 
theory of star formation in NGC 2264 (Sung et al. 11997) as de- 
scribed above. 

Another possible explanation for the differences in the two 
stars' fundamental parameters could be that VAS 20 is part of 
a binary system. We have assessed that the CoRoT photome- 
try is not affected by a contaminating source and we do not 
find any evidence for binarity in our spectroscopic measure- 
ments and in the spectral energy distribution. The log g value of 
VAS 87 derived from our spectroscopic data agrees well to the 
corresponding value determined from the evolutionary tracks. 
But for VAS 20 there is a difference of about 0.3 dex between 
the two values. We therefore hypothesize that the luminosity of 
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VAS 20 is overestimated when extracted from spectroscopy be- 
cause VAS 20 might be part of an unresolved binary system. 

Table 4. Fundamental parameters for VAS 20 and VAS 87. 





VAS 20 


VAS 87 


r eff [K] 


6380 ± 150 


6220 ± 150 


l°g£ 


4.0 ± 0.2 


3.8 ± 0.2 


v mic [kms 1 ] 


1.8 + 0.5 


1.7 ±0.5 


ysin/ [km s -1 ] 


42 ±2 


18 ± 1 


M[M ] 


2.17 ±0.15 


1.55 ±0.15 


log L [L Q ] 


1.2 ±0.1 


0.8 ±0.1 


R[R Q ] 


3.5 ± 0.5 


2.2 ± 0.5 



The low temperatures of the "comparison" stars ensure that their 
variability is likely to be produced by rotational modulation, as 
no pulsation (other than solar-like oscillations which are in a dif- 
ferent frequency range, i.e., larger than the Nyquist frequency of 
84 d~') can be present, and that their level of activity is equal 
or higher than that of VAS 20 and VAS 87. A typical example 
is given in Figure |T0j the light curve (top panel) and amplitude 
spectrum (bottom panel) of CoID0223980447 observed in 2008 
is shown. CoID0223980447 has a spectral type of K6. The fre- 
quency analysis shows two independant frequencies (Fl and F2 
in Figure ITOb and a third frequency which is 1/2 Fl indicating 
the presence of non-linear effects. From the frequency analy- 
ses of the spotted stars, we can confirm the findings by Alencar 
et al. (1201 Oi l and we infer that rotational modulation can be re- 
sponsible for a maximum of three frequencies in the 23-day long 
CoRoT data sets from 2008. 



4.2. Comparison with photometric colors 

Young stars are often still surrounded by the remnants of their 
birth clouds which affects the observed photometric colors. 
Therefore the application of photometric calibrations (developed 
for main sequence stars) to PMS objects will often lead to erro- 
neous values. 

In the present case, we used the available Stromgren pho- 
tometry (Strom et al. 1 197 11 Neri et al. 11993) and the photometric 
calibration by Moon & Dworetsky (1985) to test this hypothe- 
sis. For VAS 20 we get r eff =5800K, logg=3.6 and a metallic- 
ity fFe/H] = -0.76 and 7/ eff =6400K, logg=4.6 and a metallic- 
ity fFe/H] = -0.61 for VAS 87. Although these T eS values might 
be used as a first estimate, the logg values of 3.6 and 4.6 are, 
respectively, too low and too high for these stars and the quite 
low [Fe/H] values do not match the overall cluster metallicity of 
-0. 15 (Lynga 1987). This illustrates that the results of photomet- 
ric calibrations for the determination of the fundamental param- 
eters of PMS stars cannot be trusted and that only spectroscopic 
measurements can give reliable values. 

5. Rotation & Pulsation 

As star spots and y Doradus pulsation cause variability on the 
same time scales, i.e., in the low frequency domain, we here dis- 
cuss the nature of the observed variations in VAS 20 and VAS 87 
which we attribute to a combination of rotation and pulsation. 

Surface inhomogeneities can cause complex variability fea- 
tures in the light curves and amplitude spectra due to multiple 
spots on the surface, differential rotation or migrating star spots 
(Strassmeier 2009). Nevertheless, frequencies originating from 
rotational modulation are mostly accompanied by non-linear ef- 
fects which can be identified in the amplitude spectra as harmon- 
ics and linear combinations. 

Potential PMS y Doradus pulsators are located in the same 
region of the HR-diagram as the T Tauri stars which are known 
to show regular, semi-regular and irregular variabilities . For the 
classical T Tauri stars in NGC 2264, Alencar et al. (20101 re- 
port a mean rotational period between 3 and 4 days using data 
obtained during the same CoRoT run in 2008. 

Based on this information, we performed additional fre- 
quency analyses of cluster members with known variability 
caused by spots (i.e., the classical T Tauri stars) which have spec- 
tral types from K0 to M0. These "comparison" stars were ob- 
served with CoRoT in the same run in 2008, i.e., have the same 
data set lengths, but are much cooler than VAS 20 and VAS 87. 
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HJD - 2445545.0 [d] 
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1.0 
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Fig. 10. The "comparison" T Tauri star CoID0223980447: 
CoRoT light curve (top panel) from SRaOl (2008) and origi- 
nal amplitude spectrum (bottom panel) from to 2 d~' (bottom 
X axis), i.e., from to 23.15/jHz (top X axis), where the three 
significant frequencies are identified in black. 



Using data from the Kepler satellite, Balona et al. (120121 1 find 
three types of stars with multiple peaks in the frequency domain 
of y Doradus pulsators: stars with one or two dominant peaks 
with asymmetric light curves (their ASYM group), stars with 
two close frequencies with symmetric light curves (their SYM 
group) and stars showing multiple low frequencies (their MULT 
group). While stars in the MULT group are very likely y Doradus 
pulsators, there might be a mix of spotted stars and pulsators in 
their sample of the ASYM and SYM groups (Balona et al. 12012) . 
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From their analyses the authors also conclude that the pulsation 
periods of their stars are close to the rotational periods. 

The light curves of our two stars look rather symmetric, but 
once the peak with the highest amplitude (i.e., Fl for VAS 20 and 
F5 for VAS 87; see Figures [T]and|2]l is subtracted, multiple peaks 
at low amplitudes remain which strongly indicates the presence 
of pulsation. 

Using a v sin i value of 42 km s and a stellar radius of 3.5R 
for VAS 20, the maximum rotation period corresponds to ~4.2 d. 
We identify Fl at 0.672d~', i.e., a period of about 1.5 days, to 
be likely caused by rotation because it is the highest and most 
dominant amplitude peak in our analysis with the other frequen- 
cies lying at significantly lower amplitudes (i.e., smaller than 
1 mmag). With these parameters, the projected inclination angle 
of the rotational axis is computed to be 23 degrees. 

For VAS 87 a vsini value of 1 8km s~ with a stellar radius 
of 2.27? would correspond to a maximum rotation period of 
~6.1d. In our frequency analysis, the by far most dominant 
and highest amplitude peak in the 2011/12 data set is F5 at 
0.262d~ 1 corresponding to a period of 3.8 16 d. In the 2008 data 
the highest amplitude peak is Fl at 0.519d _1 which is about 
twice F5. As the surfaces of young stellar objects are likely 
not homogeneous, we assume that in 2008 the surface inho- 
mogeneities were not that prominent than in 2011 and that 
pulsational and rotational activities occurred on similar levels 
in 2008. In this context, we still have to remember that the 
CoRoT photometry might be affected by instrumental effects to 
a presently unknown extent and that not all the large amplitude 
differences between the data sets from 2008 and 201 1/12 might 
be interpreted astrophysically. We therefore conclude that F5 at 
0.262d~ 1 (value from the 201 1 data set) corresponding to a period 
of 3.8 16 d is likely to be caused by rotation. The corresponding 
projected inclination angle of the rotation axis is then 36 degrees 
for VAS 87. 

We conclude that rotational modulation may explain only 
few of the detected frequencies in VAS 20 and VAS 87, but the 
clear multi-periodicity of both stars shows that rotation cannot 
be the only cause of their variability and that indeed g-modes 
have been observed. 



6. Conclusions 

VAS 20 and VAS 87 were observed for about 23.4 days with the 
CoRoT satellite during the short run SRaOl. VAS 87 was re- 
observed by CoRoT during the short run SRa05 in 201 1/12 for 
about 39 days. The analyses of the high-precision photometric 
time series from space revealed ten and fourteen intrinsic fre- 
quencies between and 1.5 d -1 which we attribute to be caused 
by a combination of rotation and y Doradus type pulsation. From 
high-resolution spectroscopy we determined the atmospheric pa- 
rameters and chemical abundances for VAS 20 and VAS 87. Our 
T e s values of 6380 K and 6220 K agree, within the uncertainties, 
with those of King (11998b . 

Using our measured v sin i values of 42km s _1 and 18km s 
combined with the derived stellar radii of 3.5R Q and 2.2R for 
VAS 20 and VAS 87, respectively, we identify Fl at 0.672d _I for 
VAS 20 and F5 at 0262d- ! for VAS 87 to be likely caused by 
rotation, while the other frequencies are interpreted as g-modes 
of y Doradus pulsation. 

In this context, we also compared VAS 20 and VAS 87 to sev- 
eral additional variable stars in NGC 2264 in the spectral range 
from K0 to M0. At such cool temperatures variability can only 
be caused by spots on the stellar surface. From this analysis we 



find that rotational modulation can be responsible for three fre- 
quencies at maximum, but cannot explain the presence of ten and 
fourteen intrinsic frequencies identified in VAS 20 and VAS 87. 

The probability that VAS 20 and VAS 87 are members of 
NGC 2264 is very high due to their reported X-ray fluxes, radial 
velocity measurements and proper motions. The radial velocity 
(v T ) of the cluster NGC 2264 was derived from measurements of 
13 stars and is given as 17.68 + 2.26 km s~ by Kharchenko et al. 
(2005). From our spectra we derive v T values of 17.7 ± 1.0 and 
19.3 ± 0.9 km s" 1 for VAS 20 and VAS 87, respectively, which 
are in agreement with the cluster average. 

Our spectroscopic analyses add another crucial piece of in- 
formation: the relatively high Li abundance we measured is in 
excellent agreement with the one obtained by Sestito & Randich 
(120051 1 for other stars of similar r en - in this cluster. The Li abun- 
dances reported by King d 19981 1 are higher than the values that 
we derived but agree within the uncertainties of 0.2 dex. These 
systematic differences might occur because we used a higher r e jf 
and took into account the full structure of the line. Additionally, 
our measurements of the equivalent widths are smaller by 20 m A 
compared to King ( 1998). 

As the variability observed in VAS 20 and VAS 87 is very 
likely caused by a combination of rotational modulation and g- 
mode pulsation and as both types of variations occur on the same 
time scales, it is currently hard to unambiguously distinguish be- 
tween them. But the multiple frequencies discovered in VAS 20 
and VAS 87 provide a first strong indication for an observational 
detection of y Doradus pulsation in PMS stars which was the- 
oretically predicted (Bouabid et al. 1201 1) . These first observa- 
tional candidates for PMS y Doradus stars will be an important 
input for theory and help to continue the development of astero- 
seismic models for PMS stars. 
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